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Abstract 
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^ . We propose a simple quaternionic reformulation of Maxwell's equations for 

inhomogeneous media and use it in order to obtain new solutions in a static 
case. 
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1 Introduction 

The algebra of quaternions was applied to the study of Maxwell's equations starting 
from the work of J. C. Maxwell himself. The standard reference for the quaternionic 
reformulation of Maxwell's equations is the work [|J, where the Maxwell equations for 
vacuum were written in a simple and compact form. In this relation we should mention 
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also the earlier article [IB]. Some new integral representations for electromagnetic 
quantities based on the idea of quaternionic diagonalization of Maxwell's equations 
for homogeneous media were obtained in |J (see also []I2] , Chapter 2]). A review of 
different applications of quaternionic analysis to the Maxwell equations can be found 
in 0. In the recent work |9j with the aid of quaternionic analysis techniques the 
Maxwell equations for inhomogeneous but slowly changing media were diagonalized 
and new solutions obtained. Nevertheless even the question: how to write the Maxwell 
equations for arbitrary inhomogeneous media in a compact quaternionic form remained 
open. In the present work we propose such a reformulation and use it for obtaining 
new results in the static case. 



2 Preliminaries 

We denote by H(C) the algebra of complex quaternions (=biquaternions). The ele- 
ments of H(C) are represented in the form q = Ylk=o Qkik > where ^ G C, io is the 
unit and %\. , k — 1,2,3 are standard quaternionic imaginary units. We will use also 
the vector representation of complex quaternions: q = q$+lf > where if = Ylt=i 1 k ^ k ■ 
The vector parts of complex quaternions we identify with vectors from C 3 . The prod- 
uct of two biquaternions can be written in the following form 

V ■ Q = Po%~ <f,t> + [f Pa~t + qo~P, 

where < Jf, if > and [~p* x if] denote the usual scalar and vector products respec- 
tively. 

We will use the following notations for the operators of multiplication from the 
left-hand side and from the right-hand side 

p Mq = p ■ q and M v q = q ■ p. 
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Remark 1 The scalar product of the vectors J? and if can be represented as follows 

<J?,^>= -\{ r M + M r )~cf. 

On the set of differentiable H(C) -valued functions the Moisil-Theodoresco operator 
is defined by the expression Df = Ylk=i ^kdkf ■ In a vector form this expression can 
be written as follows 

J D/ = -div7 + grad/ + rot?, 

where the first term is the scalar part of the biquaternion Df and the last two terms 
represent its vector part. 

Let us note some simple properties of the operator D which will be used in this 
work. Let ip be a scalar function and / be an H(C) -valued function. Then 

D(<p.f)=D<p.f + <p-Df, 

and 

( o-^)f = v n(^f). (i) 



3 Quaternionic reformulation of Maxwell's equa- 
tions 

The Maxwell equations for complex amplitudes of the time-harmonic electromagnetic 
field have the form 

div(e^) = p, div(^f?) = 0, (2) 
rot if = iueE* + ~~f , rot ~~e! = —iou/ilf , (3) 
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where ~E* and H are C 3 -valued functions; lo is the frequency; e is the permittivity 
and [A is the permeability of the medium. We suppose e and \i to be two times 
differentiable complex valued functions with respect to each coordinate Xk , k = 1, 2, 3 . 
Note that they are always different from zero. 
The equations (§) can be rewritten as follows 

.. -=> grade -=> p 
div E + < , E >= - 

e e 

and 

div^+<^,^>=0. 
H 

Combining these equations with (|3]) we have the Maxwell equations in the form 

e e 

and 

DTt=<V*±,7t>+iu,e-$ + T- 

Taking into account Remark [l] we rewrite them as follows 

1 grade.—* 1 . , grade d — * p 
(D + -2 )E = — M—E - ioofiH - - 

2 s 2 & 



and 



lgrad/i.— * 1 , r gad_M _ » _4. — > 

m + _2 £ Uf = — M—H + ic^E + j 

v 2 7 2 



Due to ([J) we obtain 



and 



-^=D{y/Ji ■ 7?) + H ■ Jt = iueE* + 7 

VJ 1 
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where 



e = ■= — and J7 — 



Introducing the notations 

S = \fel^ , Tt = y/JiH and k = u^/ejl, 

we finally arrive at the following system 

(D + M^)t = -ikTt - -7= (4) 

V £ 

and 

(D + M^yyl = ikt + yfjff, (5) 

which is equivalent to @, (§). This pair of equations is the quaternionic reformulation 
of the Maxwell equations for inhomogeneous media. 

The operator D + M a , with a being a constant complex quaternion was studied 
in detail in []12 |. Note that ~t and ~j£ are constants if e and // are functions of the 



form exp(axi + bx^ + ex 3 + c?) with constant a, 6 and c. For the case when a is 
not a constant there were proposed some classes of exact solutions in ||, JlO] and JIT 



In the same articles the reader can see that the classical Dirac operator with different 
potentials is closely related to the operator D + M a . A simple matrix transform 
proposed in M turns the classical Dirac operator into the operator D + M a , where 
a contains the mass and the energy of the particle as well as the terms corresponding 
to potentials. 

In a static case ( u — ) we arrive at the equations 

{D + M^)~t — — 7=, 

and 

(D + M^)Jt = yfjff. 
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Thus we are interested in the solutions for the operator D + M , where the complex 
quaternion ~o? has the form ~ot = (grad <p)/<p and the function <p is different from 
zero. Note that due to (|l|) the operator D + M permits a complete study which can 



be found in JTJ], because it practically reduces to the operator D . In the case of the 
operator D + M the situation as we will see later on is quite different. 
Consider the equation 

(D + Af*)~? = 0. (6) 

Denote 

_ Ay? 
V> ' 

In other words (p is a solution of the Schrodinger equation 

-Atp + v(p = 0. (7) 

Proposition 2 Let ib be another solution of (Qj. Then the function 

7 = (D- (8) 

is a solution of (j^). 

Proof. The proof consists of a simple calculation. Consider 
L>7 = -Aib -Dib- — -ib- £>(— ) 

= —vip — Dtp ■ — — + ip ■ ( — —J + tp- 



<p V 3 



This proposition gives us the possibility to reduce the solution of (^j) to the Schrodinger 
equation (|7]). Moreover, if -0 is a fundamental solution of the Schrodinger operator: 
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then the function / defined by (^) is a fundamental solution of the operator D + M~^ 
that can be seen following the proof of Proposition |^. 

Remark 3 The proposition [| is closely related to the factorization of the Schrodinger 
operator proposed in fi^J, Jj]/. Namely, for a scalar function u we have the equality 

(D + M a )(D - M a )u = (-A + v)u, 

if the complex quaternionic function a satisfies the equation 

Da + a 2 = -v. (9) 

It is easy to check that for a = (grad <£>)/<£> the equation (fjj is equivalent to (ffi). 
Equation (j^) can be considered as a natural generalization of the ordinary differential 
Riccati equation. In j^J and jl^/ the corresponding generalizations of the well known 
Euler theorems for the Riccati equation were obtained. 

Let us consider the following simple example of application of the Proposition |2|. 

Example 4 Consider equation (||) in some domain ftcl 3 and let Aip/ip = —c 2 in 
Q , where c is a complex constant. In this case we are able to construct a fundamental 
solution for the operator D + M . Denote 

p ic\x\ 



47T |x| 

This is a fundamental solution of the operator — A— c 2 . Then the fundamental solution 

of D + M is constructed as follows 

-j> grady?(a;)^ e"^ . x ^ . x grady2(a:)^ e lc ' x \ 

ip(x) An \x\ |x| 2 \x\ 47r |x| ' 

where x = Yj\=i x k% ■ 

Note that it is not clear how to obtain this result for the Maxwell operators D+M~^ 
and D + using other known methods. 
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